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Recent works have suggested that selection criteria based on mid-IR properties, i.e. extreme colors and bright flux levels, can be used 
to reveal a population of dust-enshrouded, extremely-luminous quasars at z ~l-2. However, the X-ray spectral properties of these 
intriguing objects still remain largely unexplored. We have performed an X-ray study of a large sample of bright mid-IR (F 2 4,„„ > 
1.3 mJy) galaxies showing an extreme MIR/Optical flux ratio (F^^/Fs >2000) in order to confirm the presence of a luminous active 
nucleus in these very red objects. Sampling of a large area is required to pick up objects at the highest luminosities given their low 
surface density. Accordingly, we have applied our selection criteria to an area of ~6 deg 2 covered by XMM-Newton/ Chandra ob- 
servations within the ~50 deg 2 SWIRE survey, resulting in a final sample of 44 objects. The vast majority of the source redshifts, 
both spectroscopic and photometric, are in the range 0.7 i z ~ 2.5. The X-ray coverage of the sample is highly inhomogeneous 
(from snap-shot 5 ks Chandra observations to medium-deep XMM-Newton exposures of 70 ks) and, consequently, a sizable fraction 
of them («43%) remains undetected in the 0.5-10 keV band. Using spectral or hardness information we were able to estimate the 
value of the absorbing column density in 23 sources. 95% of them are consistent with being obscured by neutral gas with an intrinsic 
column density of Nh> 10 22 cirT 2 . Remarkably, we also find that ~55% of these sources can be classified as Type 2 quasars on the 
basis of their absorption properties and X-ray luminosity. Moreover, most of the X-ray undetected sources show extreme mid-IR 
colors, consistent with being luminous AGN-powered objects, suggesting they might host heavily obscured (possibly Compton-thick) 
quasars in X-rays. This demonstrates that our selection criteria applied to a wide area survey is very efficient in finding a large number 
of Type 2 quasars at z ~ 1 • The existence of this class of very powerful, obscured quasars at high z could have important implications 
in the context of the formation and cosmological evolution of accreting supermassive black holes and their host galaxies. 

C^l ' Key words. Galaxies: active - Galaxies: high-redshift - Galaxies: nuclei - Infrared: galaxies - X-ray: galaxies 
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^ 1. Introduction Intriguingly, as stressed by Fiore et al. (2003), about 20- 

30% of the sources from C hand ra/XMM-New ton surveys are 

^ ; The similarity in the anti-hierarchical evolution of both star for- optically-faint (i.e. R > 25), and with a very high X-ray-to- 

mation rate and AGN activity at z < 2 (La Franca et al. 2005; optical flux ratio (F 2 -io/F s >10). Optical spectroscopy of these 

Arnouts et al. 2007) reveals that the assembly and evolution of extreme F 2 -io/F R objects is therefore extremely challenging even 

Vj ; galaxies and those of the super massive black holes (SMBHs) with the largest telescopes. Results from deep exposures and ac- 

at their centers are intimately connected (e.g. Page et al. 2004; curat e multiwavelength studies indicate that the vast majority of 

Alexander et al. 2005; Hopkins et al. 2006; Best et al. 2005; these sources could be obscured quasars at high-z (i.e. z >1; e.g. 

Croton et al. 2006). However, the census of accreting SMBHs Mainieri et al. 2005; Maiolino et al. 2006). We have therefore 

through different cosmic epochs is still highly incomplete and started to uncover the long-sought missing AGN population of 

represents one of the major challenges in modern extragalac- the X-ray absorbed, i.e. Type 2, quasars (QS02s hereafter), that 

tic astronomy. A clear-cut evidence for the existence of a large j s a key ingredient for both the CXB synthesis models as well 

number of missing AGNs arises by the fact that the sources de- as the prediction of the SMBH/galaxy formation and evolution 

tected in the recent Chandra/XMM-Newton surveys (e.g. Brandt models (e.g. Silk & Rees 1998; Fabian 1999; Granato et al. 2004; 

& Hasinger 2005) are able to resolve only the 50-60% of the Hopkins et al. 2006). In particular, it has been hypothesized that 

cosmic X-ray background (CXB) in the 5-10 keV band (Worsley before shining in the optical as bright point-like blue objects, 

et al. 2005), and their integrated radiation significantly fails to quasars undergo an early dust enshrouded phase being associ- 

account for the CXB emission peak at ~ 30 keV (Gilli, Comastri, a ted with the rapid SMBH growth. This obscured phase is likely 

Hasinger 2007; Comastri 2004). A population of heavily ob- triggered by multiple galaxy encounters and mergers that favor 

scured AGNs showing a spectral shape consistent with the unre- the infall and the concentration in the nuclear region of large 

solved CXB component is usually invoked to fill this gap. reservoirs of gas feeding the SMBH as well as the assembly of 

the massive spheroids at high z- The link between the growth 

of the SMBHs and their host galaxies is believed to be strong 
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and bulge luminosity, mass and velocity dispersion (Merritt & 
Ferrarese 2001). A full understanding of AGN feedback pro- 
cesses is therefore a key question to be addressed in the context 
of galaxy evolution (e.g. Di Matteo, Springel & Hernquist 2005; 
Alexander et al. 2005; Menci et al. 2008; Somerville et al. 2008). 

In the latest years there have been many efforts aimed at 
unveiling the missing population of heavily obscured quasars 
at high z using alternative approaches with respect to the pure 
X-ray selection. Multiwavelength selection criteria based on 
the IR band seem to be very promising, since absorbed AGN 
optical/UV/X-ray light should be isotropically re-emitted by the 
obscuring material at these frequencies. In particular, multiple 
independent lines of evidence point towards selection criteria 
using mid-IR (3-30 pm; MIR hereafter) photometry or combi- 
nation of MIR and multiwavelenght data, to efficiently collect 
large samples of optically-faint quasars candidates at z ~ 1 (e.g. 
Houck et al. 2005; Yan et al. 2005; Weedman et al. 2006; Polletta 
et al. 2006; Alonso-Herrero et al. 2006; Daddi et al. 2007; 
Hickox et al. 2007 Donley et al. 2008 and references therein). 
Combining MIR and radio data, Martinez-Sansigre et al. (2005) 
claimed the discovery of a distant, optically-faint QS02 popula- 
tion whose number density is comparable with that observed for 
classical optically bright quasars. This might explain the paucity 
of obscured objects in the high luminosity range of the X-ray 
selected AGN population that has been reported in many works 
(Ueda et al. 2003; Steffen et al. 2003; La Franca et al. 2005; 
but see also Dwelly & Page 2006). Furthermore, recent studies 
have provided strong evidence that a large fraction of the dust- 
obscured galaxies (DOGs hereafter) with extreme F24^ m /Fj? ra- 
tios ( i 1000) and bright at 24 pm (i.e. F 2 A fl m ~ 0.5-1 mJy), 
harbor a high-z obscured active nucleus (e.g. Dey et al. 2008; 
Polletta et al. 2008a; Georgantopoluos et al. 2008; Fiore et al. 
2008a, 2008b). 

However, an accurate determination of the X-ray column 
density of the absorber in these sources is hampered by the faint- 
ness of these sources in the X-ray band (i.e. Polletta et al. 2006; 
Martinez-Sansigre et al. 2007). Therefore, most of the X-ray in- 
formation has been derived by the stacking of ChandrajXMM- 
Newton data or, for small number of sources detected with only 
a handful of counts, by the hardness-ratio analysis. 

Evidence for X-ray spectroscopically-confirmed heavily- 
absorbed QS02s at cosmological distance still remains ambigu- 
ous (i.e. Norman et al. 2002; Alexander et al. 2008; Erlund et al. 
2008). Indeed, the Compton-thick (i.e. with N H i crj 1 -1.6 x 
10 24 crrT 2 ) nature of the best Compton-thick quasar candidate, 
IRAS 09104+4109, has been recently questioned (and not con- 
firmed) by Piconcelli et al. (2007a) using XMM-Newton data. 
One of the reasons of this difficulty is that building up large 
and complete samples of bright (i.e. with a 2-10 keV flux 
F2-10 ~ 10" 14 erg cirT 2 s _1 ) QS02s by homogeneous selection 
criteria is complex and time-consuming, because of the large 
area that must be covered with deep multi-band observations. 
This is necessary due to the low space density of the AGNs 
belonging to the bright end of the X-ray luminosity function. 
The Spitzer Wide-area InfraRed Extragalactic (SWIRE) survey 
(Lonsdale et al. 2003) has imaged nearly 50 square degrees 
in 6 fields: ELAIS-S1, ELAIS-N1, ELAIS-N2, Lockman Hole, 
XMM-LSS and CDFS with medium-deep MIPS and IRAC pho- 
tometry, providing a unique opportunity to build up a large sam- 
ple of luminous QS02s. 

In this paper we report results on an X-ray spectroscopic 
study of a sample of SWIRE sources selected on the basis of 
their extreme mid-IR/optical flux ratios (F24 A , m /F« >2000) and 24 
pm flux density (i.e. F-i^m >1.3 mJy). These selection cuts have 



been applied in order to preferentially select sources harboring 
a luminous, dust-enshrouded AGN (e.g. Polletta et al. 2008a; 
Dey et al. 2008). Our aim is to confirm the presence of the ac- 
tive nucleus in these extreme DOGs and to directly measure the 
column density of the absorber via the characterization of their 
X-ray spectral properties. 

The optical magnitudes used throughout the paper are Vega 
magnitudes. The cosmological parameters used are Hq=10 km 
s" 1 Mpc" 1 , Q A = 0.7 and q = 0. 

2. The sample: selection criteria and datasets 

2.1. Optical and IR observations 

Infrared and optical photometric data are available from the pub- 
lic archive of the SWIRE survey (Lonsdale et al. 2003 and ref- 
erences therein). The catalog contains medium-deep photomet- 
ric data in MIPS 24 pm band down to 400 pJy, 5cr, and in the 
four IRAC channels from 3.6 to 8.0 pm (4.1 p}y at 3.6 pm, 5<x, 
Surace at al. 2004). Optical photometry is also available in differ- 
ent bands for selected areas in each of the 5 SWIRE fields used in 
this work (ELAIS-S1, ELAIS-N1, ELAIS-N2, Lockman Hole, 
XMM-Newton SXDS). We did not consider the CDFS here since 
these data were already studied by Fiore et al. (2008a) using a 
F24/im/FR selection technique similar to our approach. ELAIS- 
Nl and N2 have optical photometric information in 5 bands 
(U'g'r'i'Z') from the Wide Field survey (WFS, McMahon et al. 
2001). The Lockman Hole sky region is covered with 4-band 
Ug'r'i' photometry from the SWIRE photometry programme. 
XMM-LSS has very deep (R AB < 27.7) 5 band BVRi'z' photom- 
etry in 1.12 deg 2 from the Subaru XMM Deep Survey (SXDS, 
e.g. Sekiguchi et al. 2005; Furusawa et al. 2008). Full details 
on the optical and IR photometry for the ELAIS-S1 field can be 
found in Feruglio et al. (2008). 

We firstly selected sources showing a 24 pm flux density 
F(24pm) > 1000 yuJy and F 24/ jm/F« > 1000 in order to have an 
initial sample of DOGs (560 sources e.g. Fig. la. Several sources 
in ELAIS-S1, ELAIS-N1, ELAIS-N2 and Lockman Hole do not 
show an R magnitude value in the catalogs mostly since they are 
faint objects, i.e. with a magnitude above the detection limit. For 
these sources we performed aperture photometry at the position 
of the MIPS source. The photometry was performed using the 
IRAF apphot tool, with an 1.2 arcsec aperture radius (consistent 
with the released optical aperture photometry), and a 10 arcsec 
wide annulus for background estimate. For undetected sources 
we computed a 3cr limiting magnitude using the median value 
of <r rms at the positions of the source and the relation 

R Um = -2.5 x log( -y/(A) x 3 x cr rms ) + zeropoint (1) 
where A is the area corresponding to the aperture. 

2.2. Selection criteria 

We selected all the sources with available X-ray data from 
the initial sample of DOGs (i.e. F(24//m) > 1000 pJy and 
F24^m/F« > 1000 sources). The portion of sky covered by XMM- 
Newton and/or Chandra observations (for a total of ~ 6 deg 2 ) 
significantly changes across the five SWIRE fields discussed 
here, varying from -100% for the XMM-SXDS field to -5% for 
the ELAIS-N2 field. Nonetheless, since such a selection is com- 
pletely random and not based on specific target observations, it 
should not introduce any selection effect in the final sample. 

We started selecting sources with F24^ m /F« >2000 at the 
highest F24 //m flux density levels since bright 2Apm objects with 
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ID 


Name 


R 


f 3.6fim 










Log(F 24 ^ m /F s ) 


Log(^L 58/jm ) 


z 


z-type 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 



1 


J021559.46-045517.5 


25.65 


148.1 


248.8 


389.1 


636.6 


2364 


4,2 


45.14 


1.007 


SP" 


2 


J021625.91-050425.1 


24.39 


14.31 


34.45 


106.5 


400.0 


2174 


3.7 


45.96 


2.236 


PH* 


3 


J021656.59-05 1005.3 


24.31 


83.82 


146.1 


249.2 


490.0 


2083 


3.7 


45.12 


1.080 


PH* 


4 


J021734.65-051718.2 


26.36 


59.84 


111.0 


222.1 


419.7 


1322 


4.3 


45.64 


1.938 


PH* 


5 


J021743.09-044921.6 


24.20 


146 


270.2 


496.0 


871.9 


1508 


3.5 


45.13 


1.042 


PH* 


6 


J02 1749.00-052306.9 


22.36 


375.6 


714.2 


1271 


2317 


8105 


3.5 


45.67 


0.914 


SP" 


7 


J021754.45-043016.0 


23.34 


132.5 


243.9 


509.9 


1055 


2922 


3.4 


45.91 


1.780 


PH* 


8 


J021912.66-043308.5 


23.80 


27.39 


21.36 


- 


- 


1471 


3.3 


44.42 


0.754 


PH* 


9 


J021924.57-045300.0 


25.26 


30.38 


42.14 


67.16 


- 


1473 


3.9 


45.62 


1.979 


PH* 


10 


J021935. 19-050528.4 


22.70 


159.8 


142.5 


151.1 


328.6 


4095 


3.3 


44.76 


0.826 


SP* 


11 


J022003.58-045 145.6 


23.07 


123 


215.7 


380.7 


667.7 


3374 


3.4 


45.28 


1.080 


PH* 


12 


J022003.95-045220.4 


24.22 


48.1 


77.87 


177.6 


436.7 


2831 


3.7 


45.50 


1.443 


PH* 


13 


J160623.44+542301.9 


24.59 


35.85 


41.66 


43.79 


86.93 


1546 


3.6 


45.05 


1.466 


PH* 


14 


J160650.71+542510.1 


>24.30 


290.2 


258.2 


193.1 


818.8 


1768 


>3.5 


43.66 


0.22 


SP C 


15 


J160700.62+542416.8 


26.02 


18.04 


37.57 


95.66 


315.3 


1937 


4.2 


44.35 


63 +0 - 20 


PH« 


16 


J160813. 86+541942.4 


25.09 


63.38 


130.0 


270.4 


531.2 


2631 


4.0 


45.09 


1 


PH« 


17 


J160816.15+545351.1 


24.12 


64.9 


127.7 


333.1 


680.6 


2147 


3.5 


44.75 


70 +015 

"•'"-0.15 

1 


PH« 


18 


J160913.02+540509.7 


25.31 


61.84 


136.9 


292.7 


581.3 


1319 


3.8 


45.02 


PH* 


19 


J160913.28+542322.0 


23.86 


254.4 


478.6 


987.4 


1966 


4798 


3.8 


45.90 


1.400 


PH* 


20 


J160933.78+535720.2 


25.57 


50.39 


69.80 


60.08 


74.27 


1373 


3.9 


44.43 


1 


PH« 


21 


J160950.98+535836.3 


23.95 


131.9 


236.4 


406.9 


752.3 


1692 


3.3 


45.26 


1.138 


PH* 


22 


J161007. 17+540628.5 


23.34 


38.47 


90.59 


271.9 


737.3 


3284 


3.4 


45.95 


1.871 


PH* 


23 


J161017.70+535124.1 


>24.80 


17.48 


37.81 


80.36 


260.6 


1572 


>3.7 


44.37 


69+ - 20 


PH« 


24 


J161017.93+542332.1 


23.89 


36 


58.41 


135.9 


375.9 


2273 


3.5 


44.42 


63 +01 ° 

Q 77 + 0.48 

-'■ ' -0.77 


PH« 


25 


J161037.63+540603.5 


24.88 


12.05 


12.63 


- 


50.06 


1724 


3.7 


46.58 


PH« 


26 


J161051.03+535330.5 


25.39 


13.26 


20.86 


57.56 


197.4 


1560 


3.9 


44.22 


66 +015 

"•""-0.15 


PH« 


27 


J161052.74+543740.4 


23.53 


70.77 


66.55 


72.03 


136.1 


2646 


3.4 


45.21 


1.410 


PH* 


28 


J161 102.22+550551.6 


24.10 


71.53 


73.24 


121.4 


350.9 


2661 


3.6 


44.97 


0.995 


PH* 


29 


J161204.84+544149.2 


25.07 


33.74 


46.92 


42.97 


54.68 


1808 


3.8 


45.98 


*y co+0.40 

-Q-3Q 


PH« 


30 


J161306.28+543959.7 


24.91 


23.5 


54.15 


170.7 


548.9 


2922 


4.0 


44.40 


54 +!i - i!i 

"""-0.10 


PH« 


31 


J163723.94+410525.7 


23.80 


77.52 


65.69 


52.22 


133.9 


1856 


3.3 


44.94 


1.24 


PH* 


32 


J104409.95+585224.7 


23.56 


64.27 


150.4 


390.9 


1093 


4134 


3.6 


46.41 


2.54 


SP d 


33 


J104528.29+591326.6 


23.68 


32.18 


45.19 


90.94 


206.6 


2462 


3.4 


45.98 


2.31 


SP e 


34 


J104837. 82+572816.2 


23.84 


100.6 


170.3 


307.5 


625.3 


2355 


3.5 


45.53 


1.443 


PH* 


35 


J104847. 14+572337.6 


23.74 


158.1 


285.4 


451.8 


797.1 


2590 


3.5 


45.42 


1.218 


PH* 


36 


J104954.96+584429.4 


24.15 


18.45 


23.43 


58.82 


146.6 


1378 


3.3 


44.27 


76 +010 

"• '"-0.10 


PH« 


37 


J105604.84+574229.9 


22.21 


62.28 


85.63 


158.6 


438.5 


11280 


3.5 


45.59 


1.249 


PH* 


38 


J003314.61-432300.3 


24.13 


10.5 


16.69 


57.68 


142.5 


1648 


3.4 


44.46 


n 07+0.20 


PH« 


39 


J003316.92-431706.3 


24.30 


142.2 


295.1 


581.9 


1059 


2197 


3.6 


44.88 


0.689 


SP / 


40 


J003333.75-432326.9 


24.41 


48.02 


73.56 


161.1 


225.4 


1493 


3.5 


45.03 


I 24+0.20 


PH« 


41 


J003336.26-431731.7 


24.65 


18.66 


34.14 


82.01 


259.0 


2188 


3.7 


44.51 


a76 -Sl 


PH« 


42 


J003348. 18-432822.3 


>26.63 


16.17 


40.72 


98.36 


346.5 


1682 


>4.4 


44.70 


0.84!; 
1 69 +ii -* 


PH« 


43 


J0035 18.20-433414.6 


25.19 


12.06 


16.96 




142.6 


1745 


3.9 


45.38 


PH« 


44 


J003641. 47-432038.1 


23.91 


110 


224.6 


459.3 


911.7 


2746 


3.6 


45.25 


1 


PH« 



Column: (1) ID number; (2) Name from SWIRE2 catalogue; (3) Magnitude R; (4) 3.6 /im Flux density in yujy; (5) 4.5 /im Flux density in yujy; (6) 
5.8 jim Flux density in y<dy; (7) 8 jim Flux density in jiiy; (8) 24 //m Flux density in /Jy; (9) Logarithm of F 2 4 /J „,/F K flux ratio; (10) Logarithm of 5.8 
yum Luminosity in erg s -1 ; (11) Redshift; (12) Type of redshift: SP spectroscopic redshift [" this work (see section 2.3.1 for details); * Simpson et 
al. 2006; c from Sloan Digital Sky Survey Data Release 2; d Polletta et al. (2006); e Weedman et al. 2006a; / Feruglio et al. 2008], PH photometric 
redshift [ s this work (see section 2.3 for details); * Rowan-Robinson et al. 2008; ' fixed value of z\. 



large F24 /Jm /F« are predominantly MIR power-law (i.e. AGN- 
dominated) sources (e.g. Dey et al. 2008; Donley et al. 2008). We 
then gradually decreased the F24 /Jm flux density threshold until a 
statistically significant sample of objects was collected. In fact, 
fainter sources (F24^ m ~ 0. 1 mJy) selected with the same extreme 
^2AfimfpR colors, turn out to be dominated by star-formation ac- 
tivity both in the X-ray and MIR bands (Donley et al. 2008; Pope 
et al. 2008a). Accordingly we selected all the sources showing 
a F 24/ot /Fr >2000 and a 24 //m flux density F(24jum) > 1.3 mJy 



(see Fig. la) obtaining in this way a sample of 50 objects, re- 
ferred as Extreme DOGs (i.e. EDOGs) hereafter. Since in four 
cases the same optical object has a couple of IRAC-MIPS coun- 
terparts, we retained as the most likely counterpart the MIR 
source located at the smallest distance from the optical position. 

In this way we got a sample of 46 EDOGs within the 6 deg 2 
of SWIRE survey covered by X-ray observations. This implies a 
number density of 7.7 deg~ 2 that is in very good agreement with 
that obtained applying the same selection criteria to the sam- 
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Fig. 1. a) Left panel: 24 pm flux density vs. F^um/Ffi in logarithmic scales for all the sources with F24 /Jm /F« > 1000 and FiA^m > 1 
mJy in the 5 SWIRE fields considered in this work. Sources with (without) available X-ray coverage are represented with filled 
(open) circles. Dashed lines represent the selection criteria: Vi^mfpR >2000 and 24 pm flux density > 1.3 mJy. b) Right panel: An 
IRAC color-color diagram of the SWIRE Fo^,,, > 1 mJy sample (small points). Sources in our sample with (without) Nh estimate 
are plotted with open circles (asterisks). Small (large) symbols represent sources with effective exposure time < 15 ks (> 15 ks). 
The solid line marks the region of IRAC colors typical of AGN (Lacy et al. 2004). 



pie of DOGs in the «8.6 deg 2 of the NOAO Deep Wide-Field 
Survey Bootes field analyzed in Dey et al. 2008 (i.e. 7.6 sources 
deg -2 ). The fraction of sources with F24 // ,„/Fr >2000 among all 
the sources with F(24yum) > 1.3 mJy is ~10%. 

Two sources (namely SWIRE2 J160734.22+544217.3 and 
SWIRE2 J160945.43+534944.6) fall in a CCD gap in the X- 
ray image and therefore were excluded from the sample since 
there is no useful X-ray information for them. The final sample 
therefore consists of 44 sources and is listed in Table 1. Fig. lb 
shows the MIR color-color diagram (e.g. Lacy et al. 2004) for 
the sources with F24 //m >1 mJy in the 5 SWIRE fields consid- 
ered in this work (SWIRE F^fm >1 sample hereafter) obtained 
using the flux densities measured in the four Spitzer-IRAC chan- 
nelo Two distinct populations can be readily separated in this 
plot. One has blue colors in F5 8/F3 .6 and red colors in F8.0/F4.5, 
most likely consisting of a mixture of starlight-dominated galax- 
ies and star-forming galaxies with MIR 5-12/mi spectra dom- 
inated by PAH emission features at low-redshift (z < 0.2-0.4). 
The other population has redder colors and can be mainly identi- 
fied with galaxies with a MIR power-law spectrum that is a typi- 
cal signature of an AGN-powered emission. The thick lines mark 
the region of the diagram typically occupied by AGNs (Lacy et 
al. 2004; Alonso-Herrero et al. 2006; Polletta et al. 2006; Lacy 
et al. 2007; Barmby et al. 2008; see also Donley et al. (2008) for 
a detailed discussion about the level of contamination by star- 
forming galaxies in this region). It is worth noting that all but 
one sources in our sample lie within the AGN region, and most 
of them show extreme red colors i.e. log(F^^fFj,d) > 0.4 and 
log(Fg 0/F4.5) i 0.4 similarly to the luminous obscured quasars 
at z=1.3-3 in the Polletta et al. (2008a) sample. The only source 
falling outside the box has the lowest spectroscopic redshift in 



1 Source Nos. 8, 9, 25 and 43 are undetected in one (or two) IRAC 
bands and, therefore, they are not plotted in Fig. lb. 



the sample, i.e. source No. 14 at z=0.22, and therefore it is likely 
that the contribution from the host galaxy stellar light can ex- 
plain the low value of the F5.8/F3.6 ratio. The range of MIR colors 
observed for our EDOGs can be likely ascribed to the relative 
strength of the contribution of star formation to the total MIR 
luminosity (e.g. Dey et al. 2008; Sajina, Lacy & Scott 2005). 
Interestingly, Donley et al. (2007) and Cardamone et al. (2008) 
found that the most luminous X-ray sources with L2-10 i 5 x 
10 44 erg s _1 in their large samples of MIR sources fall in this 
box of the IRAC color-color diagram. 

2.3. Redshifts 

Spectroscopic redshifts are available for 7 objects in the sample 
(0.22 < z sp < 2.54; e.g. Table 1). In particular, two of them are 
obtained from new observations performed with the Subaru tele- 
scope (see Section 2.3.1 for details of data reduction and analy- 
sis). 

For the remaining 37 sources we were able to derive photo- 
metric redshifts thanks to the available optical and IR photomet- 
ric data. We used the code described in Fontana et al. (2000) and 
Fontana (2001). Photometric redshifts were calculated using a^ 2 
minimization technique on spectral libraries which include star- 
burst, passive galaxy and AGN semi-empirical template spectra 
from Polletta et al. (2007), Fiore et al. (2008a) and Pozzi at al. 
(2007). All the available photometric data points, from the opti- 
cal to the MIPS 24 pm band, were used in the computation (e.g. 
Feruglio et al. 2008). The photo-z are reported in Table 1 with 
errors corresponding to a confidence level of 1 .6 <x and their dis- 
tribution spans in the range 0.54 < z p h < 3.77. 

In four cases the x 2 minimization technique was not able to 
properly constrain the z value. For these sources (namely Nos. 
16, 18, 20 and 44) we decided to assume a fixed value of z=l 
which is close to the mean value of the z distribution in our sam- 
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Wavelength (Angstroms) Wavelength (Angstroms) 

Fig. 2. a) Left panel: Near IR spectrum of source No. 1 at z=1.007. b) Right panel: Near IR spectrum of source No. 6 at z=0.914. 
Both spectra were obtained by MOIRCS/SUBARU observations. 



pie ((z) ~1.2). It is worth noting that this is a conservative as- 
sumption since the typical value of z found for large samples 
of DOGs is ~ 2 (e.g. Polletta et al. 2008a, Dey et al. 2008 and 
references therein). 

For 19 of 37 sources photometric redshifts from the SWIRE 
Photometric Redshift Catalogue published in Rowan-Robinson 
et al. (2008) are availably Their approach is based on galaxy 
and quasar templates applied to data at 0.36-4.5 pm, and on a 
set of 4 IR emission templates fitted to IR excess data at 3.6-170 
pm. The photo-z distribution from the SWIRE Photometric 
Redshift Catalogue spans in the range 0.754 < z p h < 2.236. 

2.3.1. Spectroscopic redshifts 

We obtained the spectroscopic redshifts for two sources (namely 
source Nos. 1 and 6) via near-IR spectroscopy. These obser- 
vations were performed with MOIRCS (Multi-Object InfraRed 
Camera and Spectrograph) mounted on the Subaru telescope 
(Hawaii) on 2007 Dec 22, as backup program (program No. 
S07B0087N, P.I. M. Salvato). We used the zJ500 grism (opti- 
mized for the observation of 0.9-1.78 pm region) with a R500 
filter blocking the contamination of the second-order spectra at 
A > 1 .5 pm. Each observation consists of 2 frames dithered along 
the slit by 5 arcsec. 

After standard reduction the spectra have been subtracted 
by each other in order to improve sky subtraction and flat- 
fielding and then combined. SWIRE2 J021559.46-045517.5 has 
been observed for a total of 1200 sec and SWIRE2 J021749.00- 
052306.9 for a total of 1800 sec. Sky lines have been used for the 
wavelength calibration. Results can be summarized as follows: 

- SWIRE2 J02 15 59.46-045 5 17. 5 (No. 1): We obtained a red- 
shift of z= 1 .007 on the basis of the identification of H„/16563 
(FWHM < 1 100 km s 1 ) [NII]/16548 and [NII]^6583 narrow 
emission lines, (see Fig. 2a). 

- SWIRE2 J02 1749. 00-052306.9 (No. 6): We obtained a red- 
shift of z=0.914 on the basis of the identification of strong 
[OIII]/l4959, /15007, H ff /16563 and [NII]^6583 narrow emis- 
sion lines (see Fig. 2b). This value does not match with the 

2 For 15 sources we were able to compare the photo-z from the 
SWIRE Photometric Redshift Catalogue with those computed with our 
method. We found that for ~50%(~85%) of the sources, the photo-z 
values are consistent within 15% (30%). 



tentative redshift of z=0.987 proposed by Lacy et al. (2007), 
using very low S/N optical spectral data. The FWHM of the 
H a line is FWHM<1000 km s _1 . This support a narrow line 
quasar classification for this EDOG. 

2.4. MIR properties of EDOGs 

We computed 5.8 pm monochromatic luminosities of our 
sources using a linear interpolation between the observed 8 pm 
and 24 pm flux densities at the wavelength corresponding to 
5.8 pm at the source rest-frame (e.g. Fiore et al. 2008a for fur- 
ther details). The F24 /Jm /F« values as a function of the AL^^ m 
luminosity for the SWIRE F^Apm > 1 mJy sample and for our 
EDOG sample are compared in Fig. 3a. 

It can be seen that our selection criteria (F24 A ,,„/Fr > 2000 
and F24 /Jm > 1.3 mJy) are well suited to find highly obscured, 
high-luminosity (ALs.^m > 10 44 erg s~') sources. Furthermore, 
the combination of large F24 //m and extreme MIR colors guaran- 
tees that the vast majority of our EDOGs are AGN-dominated in 
the MIR band (e.g. Polletta et al. 2008a, Dey et al. 2008, Donley 
et al. 2008). Our finding is therefore in agreement with the corre- 
lation between F24 // „ i /Fr and MIR luminosity reported by Fiore 
et al. (2008a) for obscured AGNs (see their Figs. 2a and 2b). 

On the other hand, highly luminous (i.e. ALs.%^, > 10 44 erg 
s ) MIR sources with a F24 / ,,„/Fr value in the range 10-100 pre- 
dominantly show a continuum-dominated MIR spectrum with 
PAH emission features, suggesting they are composite ULIRG 
systems containing both a dominant AGN and a starburst com- 
ponent (Sajina et al. 2007, Yan et al. 2007). MIR-selected, lumi- 
nous Type 1 quasars typically fall in this region of the F24 jUm /Fs 
- ALs.Hf,m diagram (e.g. Brand et al. 2007; Fiore et al. 2008b; 
Polletta et al. 2008a). Finally, sources at lower luminosities (i.e. 
A Ls.Sftm < 10 44 erg s~') can be likely associated with Seyfert-like 
AGNs (obscured and unobscured) and starburst galaxies. 

In Fig. 3b is shown the distribution of the ALi^ m lumi- 
nosities as a function of the redshift for our EDOGs sam- 
ple. For comparison we also plot the sources from the 
CDFS and COSMOS fields (Fiore et al. 2008a, 2008b) with 
F24^m/Ffi >2000. The combination of our selection criteria 
with the large area of the SWIRE survey allows us to collect 
the most IR luminous sources at all redshifts: objects in our 
EDOG sample are, on average, > 0.5 ( i 1.5-2) order of magni- 
tude brighter than those selected in the 2 (<0.1) deg 2 COSMOS 
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Fig. 3. a) Left panel: AL^.^m vs. F24 /im /F« ratio for the SWIRE Fo^m > 1 mJy sample (small points). Sources in our sample with 
(without) Nh estimate are plotted with open circles (asterisks). Small (large) symbols represent sources with effective exposure time 
< 15 ks (> 15 ks). b) Right panel: The ALs.s^m - redshift plane for sources in our sample (symbol as in the left panel), COSMOS 
sources (open triangles) and CDFS sources (open squares) with F24 /Jm /FR >2000 from Fiore et al. (2008a, b). 



(CDFS) survey. Assuming that EDOGs host an AGN at their 
center, it is possible to give a rough estimate of the bolometric 
luminosity (L Bo i) using the relationship logL Bo i = logAL^^ Mm + 
1.16 from Elvis et al. (1994). In particular, Polletta et al. (2008a) 
point out that this relationship likely provides only a lower limit 
on real value of Lgoi- By proceeding in this way we obtained val- 
ues of bolometric luminosities in the range 10 12 i Lgoi ~ 10 14 
L Q . EDOGs are therefore placed at the brightest end of the bolo- 
metric luminosity function of galaxies in the Universe. 



for the MOS(PN) cameras were selected. XMM-Newton obser- 
vations sometimes experience high background flaring periods 
due to the crossing of clouds of soft (i.e. E<1 MeV) protons, 
produced by the sun and projected towards the Earth in the so- 
lar wind, that are accelerated by magnetic reconnections in the 
magnetosphere. In order to produce a cleaned event file we re- 
moved all the high background periods identified in light curves 
at energies >10 keV, where the contribution from the emission 
from real X-ray sources is minimized. 



3. X-ray observations 

3.1. Data reduction 

We searched for all Chandra and XMM-Newton archival ob- 
servations available for the SWIRE fields listed in Sect. 2.1 
as of February 2008. The X-ray coverage of the sources in 
our sample is highly inhomogeneous, ranging from shallow 
5 ks Chandra observations of the Bootes /NOAO field, up to 
medium-deep observations in the XMM-Newton SXDS field 
(~50 ks) and in the Chandra Lockman Hole region (~70 ks). 
Twenty-one sources fall within the FOV of an XMM-Newton ob- 
servation, while the sources with Chandra data are twenty- three. 

Event files of the Chandra observations were retrieved from 
the Chandra X-ray Center (CXC) via the Web ChaSeR (Chandra 
Archive Search and Retrieval) interface. The data reduction was 
performed following the standard procedures outlined in the 
Science Analysis Threads for ACIS data at the CIAO Web site0. 

Standard XMM-Newton S AS tasks epproc and emproc (S AS 
v. 7.1. 0^3 were used to linearize the PN and MOS event files. 
The event files were processed using the latest calibration files 
and cleaned up from effects of hot pixels and cosmic rays con- 
tamination. X-ray events corresponding to patterns 0-12 (Q-4-) 



3.2. X-ray Fluxes 

The count-rate in the 0.5-10 keV band (indicated as T band here- 
after) was obtained by the command sosta in the XimageQ pack- 
age. 

Using this command, we derived the number of source 
counts within a box region with a 5 (15) arcsec side for 
Chandra (XMM-Newton) data centered on the position of the 
optical counterpart of the X-ray source and also corrected them 
for the effects due to the mirror vignetting and the off-axis point 
spread function degradation. We repeated the same procedure 
to estimate the background by the extraction of the background 
counts from a source-free region around the target. We were thus 
able to infer the background-subtracted source count rate and its 
statistical significance. A 3cr upper limit is automatically cal- 
culated by sosta if the detection significance is less than 10~ 3 . 
According to this procedure, we obtained 25 positive detections 
and 19 upper limits in the T band for the 44 sources in our sam- 
ple. 

The count rate was subsequently converted into flux using 
the PIMMS0 tool. As underlying spectral model we assumed a 
power law (with a photon index fixed to F = 1 .9) modified by 
absorption from neutral gas with a rest-frame column density of 



3 http://cxc. harvard, edu/ciao 

4 http://xmm. esac. esa. int/sas/ 



5 http://heasarc. gsfc. nasa. gov/xanadu/ximage/ximage. html 

6 http://heasarc.nasa.gov/Tools/w3pimms.html 
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ID 


Exposure 


Telescope 


CRo.5-10 


Fluxo.5-10 


Log(F 2 _ 10 /F R ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


1 


18.8 


XMM-N 


25.6±4.3 


18.2+3.1 


50.1 ± 8.4 


2 


23.2 


XMM-N 


2.1±0.9 


1.8+0.8 


<2.0 


3 


29.3 


XMM-N 


<2.4 


<2.1 


<2.3 


4 


38.0 


XMM-N 


3.9±1.2 


9.4+5.1 


79.4 ± 24.4 


5 


31.8 


XMM-N 


3.5+1.8 


2.4+1.1 


0.4 ± 0.2 


6 


32.1 


XMM-N 


86.9±5.1 


111+5.5 


131.8 + 7.7 


7 


6.8 


XMM-N 


11.5+2.0 


15.6+2.7 


8.1 ± 3.4 


8 


9.0 


XMM-N 


3.2±1.6 


3.7+1.6 


3.0 ± 1.5 


9 


21.7 


XMM-N 


4.7±1.9 


6.0+2.3 


18.2 + 9.7 


10 


23.7 


XMM-N 


2.1+1.5 


2.6+1.6 


< 0.7 


11 


17.0 


XMM-N 


8.5±3.3 


6.8+2.3 


2.6 ± 1.0 


12 


17.6 


XMM-N 


17.6+3.0 


29.1+9.0 


37.2 ± 14.7 


13 


4.7 


Chandra 


<13.6 


<16.3 


<22.9 


14 


4.4 


Chandra 


<23.7 


<41.9 


<53.24 


15 


4.4 


Chandra 


<14.2 


<30.5 


<199.5 


16 


4.4 


Chandra 


<14.3 


<25.3 


<67.6 


17 


4.1 


Chandra 


<20.6 


<42.5 


<47.9 


18 


4.5 


Chandra 


<19.0 


<33.6 


<107.2 


19 


4.6 


Chandra 


10.3±4.9 


15.0+8.7 


11.7 + 5.6 


20 


4.7 


Chandra 


<23.9 


<42.3 


<173.8 


21 


4.5 


Chandra 


23.4±7.3 


30.1+6.4 


25.1+7.8 


22 


7.0 


XMM-N 


12.7+4.7 


13.4+5.2 


6.5 ± 2.4 


23 


4.4 


Chandra 


<14.2 


<29.5 


<47.59 


24 


4.8 


Chandra 


<17.9 


<37.8 


<34.7 


25 


4.6 


Chandra 


15.3+5.2 


22.1+5.5 


44.7 ±15.2 


26 


4.4 


Chandra 


<23.6 


<49.8 


<182.0 


27 


68.4 


Chandra 


<2.1 


<3.7 


<2.3 


28 


4.4 


Chandra 


<13.7 


<23.6 


<24.5 


29 


4.7 


Chandra 


<22.0 


<25.2 


<53.7 


30 


4.5 


Chandra 


<20.0 


<45.3 


<107.2 


31 


66.9 


Chandra 


1.4+0.5 


2.0+0.8 


1.5 + 0.6 


32 


67.8 


Chandra 


2.3+0.6 


3.6+0.8 


2.2 ± 0.6 


33 


66.6 


Chandra 


2.6+0.7 


4.4+0.8 


3.0 ± 0.8 


34 


4.7 


Chandra 


6.3+3.7 


10.4+6.1 


<8.5 


35 


4.6 


Chandra 


<14.4 


<25.5 


<19.5 


36 


14.2 


XMM-N 


9.1+4.1 


5.2+1.3 


5.2 ± 2.4 


37 


4.7 


Chandra 


<14.2 


<31.9 


<6.3 


38 


17.6 


XMM-N 


<3.0 


<4.3 


<4.8 


39 


21.0 


XMM-N 


37.0+4.8 


32.5+9.7 


41.7 + 5.4 


40 


21.1 


XMM-N 


30.6+4.3 


20.0+2.1 


14.1+2.0 


41 


26.8 


XMM-N 


21.9+3.4 


25.1+5.4 


45.7 + 7.1 


42 


14.5 


XMM-N 


<2.0 


<2.7 


< 2.5 


43 


14.1 


XMM-N 


46.3+6.8 


20.4+3.0 


38.0 + 5.6 


44 


22.0 


XMM-N 


49.7+5.4 


26.5+4.2 


19.1 ±2.1 



Column: (1) ID ; (2) Total effective exposure in ks; (3) Telescope (XMM-Newton or Chandra ); (4) Background-subtracted count rate in the T 
band in units of 10~ 4 Count/s; (5) Flux in the T band in units of 10~ 15 erg s -1 cirT 2 ; (6) Logarithm of the F2-10/FR ratio. 



N H =10 23 cm 2 . This is an intermediate N H value for F 2 4 //m /F s 
selected DOGs (see below). 



Table 2 lists the resulting T band count rate and flux along 
with the effective exposure time for each source in the sample. 
The latter has been obtained from the nominal exposure time 
of the observation after correction for the vignetting factor at 
the source position on the detector. We found that ~45% of the 
sources have an effective exposure time that is >15 ks. Fluxes in 
the T band span froml.8xl0" 15 tol.l xl0~ 13 erg s 1 cm 2 with 
a distribution mean value of ~ 2 xl0~ 14 erg cirT 2 s _1 . 



4. X-ray spectral properties 



In order to characterize their spectral properties (with particu- 
lar emphasis on the absorbing column density N H ), we selected 
the sources showing a detection in at least one of the follow- 
ing energy bands: soft (i.e. 0.5-2 keV, S band hereafter) or hard 
(i.e. 2-10 keV, H band). This criterion indeed allows to derive 
the hardness ratio (HR) and, therefore, provide an estimate of 
the absorption column density even if just a handful of counts 
are collected during an observation. This condition was satisfied 
by 23 out of 44 sources, i.e. ~50% of the sample and, for the 12 
brightest sources, i.e. with >40 net counts in the T band, we were 



8 



G. Lanzuisi et al.: Revealing obscured Quasars in the SWIRE survey 




10 s - 



10 23 



10 22 



10 21 - 



In 



QS02 



10 22 10 2! 
N H from HR (en 



10" 



10 4 ' 



10 4 



^-io^v (erg s-') 



Fig. 4. faj Le/f panel: Comparison of the best fit values of Nh to those obtained with HR technique, for the 12 sources with 
available spectral analysis, (b) Right panel: Distribution of column density Nh as a function of the Lq-w. Filled squares indicate 
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Nh parameter. Sources with a detection only in the Hard band have a lower limit both in the Nh and L2-10 values. Circled points 
represent the spectroscopically-identified sources. Dashed line marks the locus for sources with Nh and L2-10 values typical of 
QSQ2s (N H > 10 22 cm ~ 2 , L 2 _ 10 > 10 44 erg s" 1 ). 



also able to perform a basic spectral analysis (we refer hereafter 
to this sample of 23 EDOGs as the 'X-ray sample'). 

4. 1 . Spectral Analysis 

All sources detected with a number of T band counts >40 are 
observed with XMM-Newton. Spectra were acquired from the 
cleaned event files of each EPIC camera, i.e. PN, MOS1 and 
MOS2. 

Source counts were extracted from circular regions with a ra- 
dius in the range of 10-20 arcsec. The radius was chosen in order 
to optimize the signal-to-noise ratio (S/N) for each source, but 
in same cases it was limited by the presence of nearby sources 
or CCD gaps. Appropriate response and ancillary files for all 
the EPIC cameras were created using RMFGEN and ARFGEN 
tasks in the XMM-Newton SAS, respectively. In case of multi- 
ple observations of the same source, we combined the source 
and the background spectra from different exposures using the 
tool mathpha in order to improve the photon statistics. In doing 
so we also computed the averages of Response Matrix Files and 
Ancillary Region Files, weighted for exposure times, using the 
tools addrmf and addarf. We performed the spectral analysis for 
12 sources having net T band counts in the range 40-600 (with a 
median value of ~80 counts). All the fits discussed in this paper 
include the absorption due to the line-of-sight Galactic column 
density (Dickey & Lockman 1990). 

For the nine sources with a number of counts <150 we 
used the C-statistic minimization method (Cash 1979), being 
particularly well-suited for low-count spectra. This is a maxi- 
mum likelihood method that allows the use of unbinned data. 
Since it assumes that the counts collected in a given channel fol- 
low a Poisson distribution, it cannot be applied to background- 
subtracted data. Therefore, before performing the spectral mod- 
eling we accurately parametrize the background spectrum, ex- 



tracted from a large source-free region for each source, with a 
model consisting of a broken power-law plus Gaussian lines, as 
suggested by Lumb et al. (2002). Then we included this back- 
ground model in the source spectral fit (fixing the values of the 
parameters and rescaling the normalizations to the area of the 
source extraction region). We assumed an absorbed power-law 
as spectral fitting model. Due to the low number of counts we 
decided to fix the power-law photon index Y to 1.9, i.e. the av- 
erage value found for large samples of quasars in the 2-10 keV 
band (e.g. Reeves & Turner 2000; Piconcelli et al. 2005 and ref- 
erences therein) 

For the remaining three sources (namely source Nos. 6, 39, 
44) with > 150 net spectral counts we were able to perform the 
spectral analysis using the x 2 minimization method. In partic- 
ular, the large number of collected counts for the source No. 
6 (SWIRE2 J021749.00-052306.9) allowed to significantly de- 
tect an Fe Ka emission line at ~6.4 keV in this high-luminosity 
(L2_io~ 5 x 10 44 erg s -1 ) obscured quasar at z ~ 1 (see Sect. 6.2 
for further details). Table 3 reports X-ray counts and spectral 
parameters for these 12 sources (labeled with SP). We also re- 
port the observed L2-10, the absorption corrected L2-10 an d the 
ratio /IL5 fym/Z^-io value for each of them. 

4.2. Hardness Ratio Analysis 

For the 11 sources detected with <40 T band counts we cal- 
culated the background-subtracted count-rate in the S and H 
band. We computed the HR value, defined by the formula C#- 
Cs/C/z+Cs, where Cs (Ch) is the net count-rate in S (H) band. 
The HR provides an indication about the flatness of an X-ray 
spectrum and, therefore, it can be used for an estimate of the 
strength of X-ray obscuration, since the photoelectric absorption 
mainly affects the soft X-ray portion of the spectrum (i.e. photo- 
electric absorption cross-section cr(E) oc E~ 3 ). 
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Table 3. X-ray spectral properties of the 23 sources showing a detection in at least one of the two energy bands, S and H. 
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ID CRo.s-2 CR 2 _io HR Flux 2 _ 10 N„ log(L x , o4s ) log(L x ) log(iL 5 . 8 /L 2 _ 10 ) Analysis 

(1) (2) (3) (4) (5) (6) 

1 19.4+3.4 6.1 +2.5 -0.52+0.19 l.ltht 

4 1.7 +0.8 2.2 +0.9 0.13+0.31 8.5^ 19.9!™ 

5 2.7 +1.3 1.8 +1.2 -0.25+0.53 2.2+^ <14.3 

6 21.1 +2.5 63.7+4.3 0.49+0.07 H3.rJ™ 10.8+|| 

7 4.0+2.8 7.4+3.9 0.30+0.44 14.2™ 25.0+g-g 

8 <6.9 4.8+4.1 >-0.53 3.4 + J| >0.1 

9 1.8 +1.6 2.8 +1.9 0.21+0.54 5.3«| 25.0+«°° 

10 1.0+0.9 <1.5 <0.25 <2.3 + ]j <6.6 



11 5.0 +2.3 3.5 +2.4 -0.18+0.39 5.8™ 4.9!»' 

12 1.9 +1.6 5.7 +2.5 0.50+0.44 27.H| 267.9 + _f^ g 
19 2.2+1.3 8.6+4.4 0.60+0.62 12.6+|| 25.0™ 

21 13.0+5.3 10.6+4.9 -0.10+0.31 24.8™ 4.0™' 

22 5.9 +3.0 6.7 +3.7 0.06+0.37 13.0*ff° 
25 10.9 ±3.7 4.5 ±5.0 -0.42+0.43 18.6™ <45.0 

31 <2.2 0.9 ±0.4 >-0.15 1.7*$ >0.1 

32 0.7 ±0.3 1.8 ±0.5 0.49±0.31 3.1^ 75.0™ 

33 0.7 ±0.3 1.9 ±0.6 0.48±0.28 3.8+°* 60.0™ 
36 3.7 ±2.5 9.4 ±3.8 0.43±0.40 4.3™ <28.0 

39 11.2 ±2.8 25.9 ±3.9 0.40±0.14 29.9 + J 3 \ 6.T\i 

40 23.2 ±3.5 6.9 ±2.3 -0.53+0.16 9.3+fj' <1.5 

41 4.6 ±1.9 17.5+2.8 0.58±0.18 23.7~™ 

43 40.7 ±5.8 3.3 ±3.0 -0.81+0.18 12.0™ <0.3 

44 31.2 ±4.1 18.4 ±3.6 -0.26±0.11 19.5+^ 1.7^ 



(7) 


(8) 


(9) 


(10) 


43.68 


43.72 


1.4 1 


SP 


43.90 


44.31 


1.31 


SP 


42.93 


43.03 


2.1 


SP 


44.49 


44.71 


0.95 


SP 


44.40 


44.48 


1.43 


HR 


43.13 


>43.13 fl 


1.29 


HR 


44.08 


44.15 


1.47 


HR 


42.82 


42.87 


1.76 


HR 


43.40 


43.57 


1.73 


SP 


43.40 


45.16 


0.34 


SP 


44.10 


44.92 


1.76 


HR 


44.19 


44.21 


1.05 


HR 


44.35 


44.39 


1.56 


HR 


45.26 


45.30 


1.33 


HR 


43.24 


>43.24 fl 


1.70 


HR 


44.10 


44.53 


2.23 


HR 


44.10 


44.59 


1.82 


HR 


43.01 


43.06 


1.21 


SP 


43.65 


43.85 


1.03 


SP 


43.88 


43.90 


1.34 


SP 


43.88 


44.05 


0.46 


SP 


44.44 


44.45 


1.04 


SP 


43.94 


44.01 


1.24 


SP 



Column: (1) ID number; (2) Count rate in the S band in units of 10~ 4 count/s; (3) Count rate in the H band in units of 10~ 4 count/s; (4) Hardness 
Ratio between the H band and the S band computed with the formula H-S/H+S; (5) H band flux in units of 10~ 15 erg s cm 2 ; (6) Rest-frame 
column density N H in units of 10 22 cirr 2 ; (7) Logarithm of the observed luminosity in the H band (erg s _1 ); (8) Logarithm of the absorption- 
corrected luminosity in the H band (erg s" 1 ); (9) Logarithm of the AL S .s^m/^-io ratio; (10) Analysis method: SP for spectral fitting, HR for 
hardness ratio. For Chandra data the 2-10 keV count rate and flux was extrapolated from the 2-8 keV value assuming T=1.9. 
" For these sources we computed the absorption-corrected luminosity using the lower limit on the Nh value. Accordingly, the intrinsic luminosity 
shown here is also a lower limit. 



To evaluate the absorbing column density Nh , we created a 
simulated spectrum adopting as model a power law (with photon 
index F = 1.9) modified by a rest-frame absorption. We varied 
the column density Nh in the range 20<log(NH)<24 cirr 2 by 
steps of Alog(Nn)=0.2, and the redshift in the range 0<z<3 by 
steps of Az=0.5. By computing the HR of the simulated spec- 
trum for each value of Nh and z, we created a set of curves cor- 
responding to different values of Nh in the plane z versus HR that 
allowed us to evaluate the absorption Nh for each source. Thus, 
we were able to compute both the observed and the absorption 
corrected Z^-io- For 2 source (Nos. 8 and 31) with a positive de- 
tection only in the H band, we obtained a lower limit on HR and, 
hence, Nh- The latter was used in the computation of the ab- 
sorption corrected luminosity Table 3 reports X-ray counts and 
spectral parameters for these 1 1 source (labeled HR). 

Previous works (e.g. Perola et al. 2004) suggested that the 
Nh values obtained through the HR technique tend to be over- 
estimated. In order to check the reliability of the HR technique 
in case of a small number of counts, we compared the Nh in- 
ferred by spectral analysis with the Nh value obtained using 
the method outlined above for the subsample of 12 sources (see 
Sect. 14.11 ). The values of Nh determined for each sources are 
in good agreement within errors as shown in Fig. 4a. Notably, 
the only source whose Nh estimate does not match the one-to- 
one relationship is the Compton-thick QS02 candidate S WIRE2 
J022003.95-045220.4 (i.e. source No. 12, see Appendix for fur- 



ther details). This suggests that, due to the intrinsically different 
spectral shape of Compton-thick sources, the HR analysis can 
lead to underestimate the Nh value in these sources. 

4.3. Absorption properties 

All the sources in the X-ray sample show a hard X-ray luminos- 
ity value typical of AGNs, i.e. L,2-\o~s> 10 42 erg s . We found 
that the vast majority (i.e.~95%) of the 23 EDOGs with X-ray 
spectral data are consistent with being obscured by neutral gas 
with an intrinsic column density of Nh> 10 22 cirr 2 (Fig. 4b). 
Unfortunately, the combination of large z values and, mostly, 
the limited photon statistics do not allow us to test the pres- 
ence of warmly ionized absorbing gas in our sources. In partic- 
ular the fraction of X-ray obscured EDOGs showing a cold ab- 
sorption column density strictly larger than 10 22 cirr 2 is ~55%. 
The uncertainties on this fraction is due to the large errors on 
the Nh value for the sources detected with a small number of 
counts. Table|4]lists the source breakdown according to the level 
of absorption. This table, as well as Figure 4b, also reveals the 
common presence (~55%) of QS02s among the EDOGs in the 
X-ray sample. This finding highlights the efficiency of our MIR 
selection criteria in collecting a large number of these elusive 
AGNs, which are easily missed in optical and X-ray surveys. 
This makes promising their application to large samples of X-ray 
well-exposed EDOGs with the goal of studying the properties of 
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Fig. 5. a) Left panel: F2-10/FR ratio as a function of the 2-10 keV luminosity for the X-ray sample. The solid line is the best 
log(F2-io/FR)-log(L2-ioi:eiO linear regression obtained by Fiore et al. (2003) for a large sample of optically obscured AGNs (see 
Sect. 4 for details), b) Right panel: F2_io/F« ratio versus column density N H for the sources in the X-ray sample (large points). For 
comparison we also plot the values obtained for a sample of the 99 brightest, hard X-ray selected sources in the wide-area (6 deg 2 ) 
XMM-LSS survey from Garcet et al. (2007) (small circles). 



the absorption at the brightest end of the AGN X-ray luminos- 
ity distribution. It is also worth noting that the sources with the 
highest absorption column densities fall within the QS02 locus 
in the Nh versus 2-10 keV luminosity plane. 

To further test the real capability of our MIR selection cri- 
teria in gathering QS02s, we compared these results with those 
obtained for a control sample consisting of 20 randomly-selected 
sources from the same SWIRE fields taken into account in our 
work. We picked out these sources for having a hard X-ray flux 
range (F2-io~10~ 15 -10~ 13 erg crrT 2 s~') and a redshift distribu- 
tion similar to those of our X-ray sample, regardless of their MIR 
and optical properties. Photometric redshifts for these sources 
are taken from Rowan-Robinson et al. (2008) (13 sources) or 
have been derived from optical and IR photometric data (7 
sources). Spectral parameters were derived by the individual X- 
ray spectral analysis of each source. The fraction of sources in 
the control sample with Nh > 10 22 cm~ 2 is about 15%. Finally, 
the fraction of QS02s is -5%. 

4.4. X-ray to Optical flux ratio 

Optically-selected AGNs typically show an F2-10/FR flux ratio, 
in the range 0.1-10, while ~ 20% of AGNs selected in hard 
X-ray deep surveys has an F2-10/FR > 10 (see Alexander et 
al. 2001; Fiore et al. 2003). Follow-up studies of these EXOs 
found that most of them are dust-enshrouded AGN at z ~ 1 
(e.g. Maiolino et al. 2006 and references therein). Fiore et al. 
(2003) discovered the existence of a striking correlation be- 
tween F2_io/F s an d 2-10 keV luminosity for non Type 1 AGNs 
whereby the most luminous X-ray sources tend to show the 
highest F2-10/FR values. This is due to the fact that Compton-thin 
AGNs have their nuclear optical/UV light strongly attenuated 
differently from the X-ray light and, therefore, their F2-io/Fs 
represents the ratio between the X-ray flux from the AGN and 
the host galaxy star light, which spans a small luminosity range. 



Table 4. The source breakdown according to the absorption 
properties. See Sect. 4.3 for details. 



Source 


X-ray 


Control 


Type 


Sample 


Sample 




(1) 


(2) 


N H < 10 22 


13-52% 


85% 


N H > 10 22 


48-87% 


15% 


QS02s 


43-60% 


5% 


N. of Sources 


23 


20 



Column: (1) Percentage of objects for each source type found in our 
sample (i.e. for the absorbed sources we report the maximum and 
minimum value derived assuming the best fit value and the lowest 
value at 90% confidence level for the N H parameter, respectively); (2) 
Percentage of objects for each source type in the control sample (20 
objects). 



Furthermore, the F2-10/FR value is a proxy for the absorption in 
AGN (i.e. the higher is F2-10/FR, the higher is Nh) and, therefore, 
can be used to select obscured AGNs. 

In order to investigate the effects of our MIR/Optical se- 
lection on the distribution of Ft-w/Fr, we computed this ra- 
tio for the sources in our sample^ (see Table 2). In Fig. [5^ we 
plot the distribution of F2-10/FR as a function of the 2-10 keV 
luminosity for the X-ray sample. A clear correlation between 
F2-10/FR and hard X-ray luminosity is found, even if character- 

7 For the sources which are undetected in the 2-10 keV band we de- 
rived an upper limit in the F2-10/F/; value extrapolating the H band flux 
using as underlying spectral model an absorbed (Nh=10 23 cm~ 2 ) power 
law with T=1.9, e.g. Sect. 3.2. 
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ized by a quite large scatter, in particular for the most X-ray 
absorbed sources. The solid line in the plot represents the best 
log(F2-io/FR)-log(L2-ioA: CT ) linear regression obtained by Fiore 
et al. (2003) from a large sample of X-ray selected, optically- 
obscured AGNs. 

The F2-io/F« value as a function of the column density 
Nh for our X-ray sample is displayed in Fig.[5t>. Notably, all but 
two sources with F2-io/F# >10 show X-ray absorption column 
densities consistent with being Nh i 10 22 cirT 2 . We also report 
in this diagram the values obtained for a sample of the brightest 
99 hard X-ray sources detected in the XMM-LSS survey (e.g. 
Garcet et al. 2007). This survey covers ~ 6 deg 2 and, therefore, 
is one of the surveys with the widest coverage area in X-rays. 
This sample of 99 optically identified sources (half of sources 
show az > 0.5) covers roughly the same range of 2-10 keV 
flux (4xl0~ 15 < F 2 -\o < 2 x 10~ 13 ergs cirT 2 s~') of our X-ray 
sample, and the Nh values were obtained from individual spec- 
tral analysis. From this figure it can be seen that EDOGs have 
higher values of F2-io/Fj? and that also the Nh values are signifi- 
cantly higher, as expected on the basis of Fiore et al. (2003) and 
Cocchia et al. (2007) results. 

5. Discussion 

5.1. Efficiency of our MIR selection criteria to collect QS02s 

Our study was aimed at testing the capability of selection cri- 
terion (F24 /Um /F s >2000 and F(24yura) > 1.3 mJy) to isolate lu- 
minous, X-ray absorbed quasars at z ~ 1 . The discovery of a 
QS02 in ~55% of EDOGs in the X-ray sample is the most in- 
teresting result of our work (e.g. Table 4). It is worth stressing 
that all the EDOGs with column densities consistent with being 
Compton-thick fall in the QSO regime, i.e. with a L2-io> 10 44 
erg s _1 . 

We have thus tested the efficiency of our MIR technique to 
recover a large number of these elusive objects missed in tra- 
ditional surveys. They remained poorly-sampled even by blind- 
searches performed in the deepest X-ray surveys (e.g. Fabian 
2002; Tozzi et al. 2006). Our choice of sampling at F(24//m) 
i 1 mJy appears crucial to identify very luminous AGN- 
powered objects (as confirmed by the L2-10 >10 43 erg s mea- 
sured for the sources in the X-ray sample) and minimize any 
strong contamination by non-AGN emission (e.g. Weedman et 
al. 2006; Dey et al. 2008, Sacchi et al. in prep.). 

Our findings also highlights the excellent opportunity rep- 
resented by this simple approach in enabling many important 
studies concerning the nature of this dust-enshrouded, bright 
AGN population at high redshifts. In particular, the systematic 
application of this method to wide-field MIR surveys with ap- 
propriately deep and homogeneous X-ray coverage will provide 
useful insights on largely unexplored issues such as: (i) the de- 
mography and surface density of highly obscured quasars; (ii) 
the multiwavelength SED of distant QS02s; (iii) the Nh distri- 
bution at the largest X-ray luminosities; (iv) the luminosity func- 
tion of Compton-thin and Compton-thick quasars at z > 1-2, i.e. 
near the peak of the AGN activity. Our work may be, therefore, 
a pathfinder study for the next-generation observatories, when 
large areas of the sky will be surveyed at faint hard X-ray fluxes. 

Our study represents the first attempt to characterize in 
details the X-ray spectral properties of EDOGs. As expected 
on the basis of previous works based on analogous (although 
less severe) MIR selection criteria (Fiore et al. 2008a, 2008b; 
Georgantopoulos et al. 2008; Polletta et al. 2008a) these ex- 
tremely red sources are found to harbor X-ray Type 2 AGNs. 



However, this MIR bright galaxy population deserves further 
and deep investigations in X-rays since it constitutes by far one 
of the most intriguing recently-discovered class of extragalactic 
objects. Extraordinary bolometric luminosities, high dust con- 
tent and redshift distribution of these sources lead to believe that 
the EDOG phenomenon can be potentially associated with the 
early stages of quasars evolution as envisaged by Silk & Rees 
(1998). 

Interestingly, as stressed by Dey et al. (2008), the space den- 
sity of the F(24 pm) > 1 mJy and ¥ 1Alim [¥ R >1000 DOGs is com- 
parable to that of unobscured quasars derived from the MIR 
luminosity function at z=2 (Brown et al. 2006). The brightest 
DOGs likely are AGN-powered sources deeply obscured both 
in the optical and soft X-rays and, therefore, cannot be revealed 
by traditional optical/soft-X surveys, which miss a sizable frac- 
tion of the accretion history of the Universe. 

In this paper we focus on the most extreme DOGs (represent- 
ing -10% of the ¥{2Apm)>\ mJy and F 24/Jm /F R >1000 DOGs) as 
a first step in the exploration of the X-ray spectral properties of 
this class of MIR bright, dust enshrouded galaxies. We indeed 
found that the vast majority of the EDOGs harbor an obscured, 
intrinsically luminous X-ray AGN. It would be very interesting 
to extend this study to the whole DOG population, in order to 
provide an accurate estimate of the AGN space density at high z. 

5.2. MIR luminosity versus Hard X-ray luminosity 

In Fig. 6a we compare the 2-10 keV luminosity with the 5.8 
pm luminosity for sources at the high-luminosity end of the 
AGN population. Objects from our X-ray sample are indicated 
with filled circles, while the shaded area represents the region 
of L2-\o-AL 5 $,f lm plain derived for AGNs in the local Universe 
by Lutz et al. (2004). In the figure are also plotted data taken 
from other samples of QS02s available in the literature. In par- 
ticular, we have included: (i) three QS02s from the Polletta 
et al. (2008a) sample of z >1 objects with F24 > 1 mJy, 
F 2 WF« >400, A L 5 . %im > 10 12 L and 10 21 < N H < 10 23 
cirT 2 (namely LH_02, LH_A5 and LFLA6; filled squares), (ii) 
six Compton-thick z ~ 2 QS02 candidates from Alexander et al. 
(2008) which are extremely weak (or even undetected) in deep 
X-ray exposures but show evidence for strong AGN activity in 
the optical and/or MIR band (filled pentagons), and (iii) eight X- 
ray selected QS02s (Zo.5-10 >10 44 erg s" 1 and N H >10 22 cm" 2 ) 
from the Sturm et al. (2006) sample (filled triangles). 

We have also included in the plot two well-known, very 
powerful (L/k > 10 12 L ) dust-enshrouded QS02s, i.e. IRAS 
09104+4109 (empty star, z=0.442; e.g. Hines & Wills 1993; 
Piconcelli et al. 2007a) and 4C 23.5ffl {empty triangle, z=2.48; 
e.g. Johnson et al. 2006; Seymour et al. 2008). Sources ab- 
sorbed by a column density of Log(Nn)>24 cirT 2 are marked 
with empty squares. 

The solid line in this plot represents the linear regression de- 
rived for our sample: 

Log{Li- W ) = 0.805 x Log(A L 5 . 8/im ) + 7.785 (2) 

once 3 sources with L2-io< 10 43 erg s _1 were excluded, since 
their weakness in the hard X-ray band (see Table No. 3) sug- 

8 We re-analyzed the XMM-Newton data of 4C 23.56. Our best-fit (an 
absorbed power law with N H = 4.4+" x 10 23 cnT 2 and T fixed to 1 .9 plus 
a soft (r =3) unabsorbed power law) is consistent with those reported 
by Johnson et al. (2006), confirming the QS02 nature of this high-z 
radiogalaxy. We derived a 2-10 keV intrinsic(observed) luminosity of 
L 2 -io=1.3(0.4)xl0 45 ergs- 1 . 
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Fig. 6. a) Left panel: ALs.^m versus L2-10 for sources in our X-ray sample (filled circles), from Polletta et al. (2008a, filled squares), 
Alexander et al. (2008, filled stars), Sturm et al. (2006, filled triangles). Empty star represents the QS02 IRAS 09104+4109 
(Piconcelli et al. 2007a) and the empty triangle represents the QS02 4C +24.36 (Johnson et al. 2007). See Sect. 5.2 for further 
details. Empty squares indicate sources absorbed by a column density of Log(Ne)>24 cirT 2 . Solid line shows the linear regres- 
sion obtained for sources with L2-10 > 10 43 erg s _1 in our X-ray sample. Short dashed line is the linear regression obtained 
from COSMOS Compton-thin AGNs (Fiore et al. 2008b). The long dashed line is the combination of the relationships AL(5100 
A)/AL(6Jfim) and L2-io//l£(5100 A) from Maiolino et al. (2007). b) Right panel: ALs.^ m versus the observed L2-io- Downward- 
pointing arrows indicate the 3<x upper limit on ^2-10 ^ or tne X-ray undetected sources (large symbols indicate sources with exposures 
>15 ks, see Sect. 6.2). Other symbols are as in the previous panel. The solid line represents the linear regression obtained from our 
X-ray sample, while the dotted and dashed lines represent the ratio L2-io/ALs.s fnn expected for an absorption of Nh=10 23 and 10 24 
cm -2 , respectively. In both panels the shaded area indicates the range of L2-m-AL^,% flm found for AGNs in the local Universe (Lutz 
etal.2004). 



gests that their ALs.^m is starburst-dominated. The long dashed 
line represents the result (i.e. Log(L2-io)=0.88xLog( AL(,j^ m ) + 
4.42 ) obtained combining the relation between the rest-frame 
continuum luminosity AL 5100 ^ and the MIR luminosity AL^j^n 
with the relation between L2-10 and AL 5m ^, both taken from 
Maiolino et al. (2007) for a sample of 25 high luminosity (i.e. 

^5iooA > er § s ~') q uasars at 2 < z < 3.5. The L2-10- 
ALs.^m relationship derived by Fiore et al. (2008b) for a large 
sample of Compton-thin AGNs in C-COSMOS and CDFS sur- 
veys (i.e. log(L2-\o) = 0.72xlog(AL^^ m )+l 1.75) is marked with 
the short dashed line. 

As expected, we found a correlation between the two lu- 
minosities. In particular this relationship is not linear, and the 
AL^.^mlLi-u) ratio increases with increasing XLs.% im . The shift 
at the highest MIR luminosities from the local Universe relation- 
ship suggests a possible evolution of the AL^ ^m/Li-w with red- 
shift. Nonetheless it is evident from this plot that there is a sig- 
nificant scatter in this relationship. This probably means that we 
have put together a mix of different classes of sources, charac- 
terized by different properties and dominant energy mechanisms. 
The line resulting from Eq.|2]divides the AL5 .^m-La-io plane in 
two regions. The region above this line includes all the X-ray 
selected QS02s from Sturm et al. (2006) and the locus typical 
of AGNs in the local Universe, which was derived by Lutz et al. 
(2004) from a sample of optically/X-ray bright Seyfert 2 galax- 



ies. Selection effects therefore favor this region to be inhabited 
by the X-ray/optically 'loud' objects. All the sources from the 
Alexander et al. (2008) fall in this part of the diagram reflecting 
the fact that they are optically bright enough to be robustly iden- 
tified via optical spectroscopy even if they likely are Compton- 
thick QS02s. The intrinsic X-ray luminosities of these sources 
was indeed estimated on the basis of the luminosity of optical 
emission lines. Two out of three objects taken from the Polletta 
et al. (2008a) sample showing moderate absorption in X-rays 
(Nh ~ 10 23 cirT 1 ) are located in this area. Interestingly, Sturm 
et al. (2006) emphasize the absence of dust emission and absorp- 
tion features in the MIR spectra of their QS02s unlike ULIRGs, 
which are very dusty systems with powerful circumnuclear star- 
bursts. This points out that these highly luminous X-ray quasars 
were able to sweep away dust envelopes from their nuclear en- 
vironment. 

In the context of the evolutionary sequence described in 
the previous Section, it seems possible that most of the objects 
falling in this region of the L,2-w-ALs l % tim plane have already 
passed through the initial dust-enshrouded phase of quasar ac- 
tivity/galaxy assembly (e.g. Silk & Rees 1998) and are going to 
evolve into classical 'blue' QSOs traditionally collected in opti- 
cal/UV surveys. IRAS 09104+4109 and 4C +24.36, which are 
respectively IR- and radio-selected, lies outside the locus of lo- 
cal Universe AGNs. 
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Table 5. The source breakdown for the X-ray sample, the Red 
sample, and the total sample. See Sect. 5.2 for details. 



Source 
Type 


X-ray 
Sample* 
(1) 


Red Sample 

(2) 


Total 
Sample 

(3) 


N H > 10 22 


50% 


93% 


66% 


CT AGNs 


13% 


27% 


18% 


QS02s 


43% 


27% 


37% 


N. of Sources 


23 


15 


38 



Percentage of objects for each source type found in the: (column 1) 
X-ray sample (e.g. Sect. 5); (column 2) the X-ray undetected Red sam- 
ple (i.e. sources showing log(F5.8/F3.6) > 0.2 and log(Fg/F45) > 0.2); 
(column 3) the combined (i.e. X-ray+undetected Red) sample. 'For the 
X-ray sample we considered the conservative values derived using the 
lowest value in 90% confidence level for the N H parameter. 



On the contrary, the majority of the sources located below 
the solid line in Fig. 6a might be objects with dust-absorbed 
powerful (ULIRG-like) nuclei in which the star-formation 
activity could significantly contribute to their MIR emission 
or characterized by intrinsic extreme MIR/X luminosity ratios 
(e.g. Polletta et al. 2008b; Piconcelli et al. 2009 in preparation). 
Accordingly, these sources may be the tip of a larger population 
representing the link between X-ray-selected QS02s with 
featureless MIR spectra (e.g. Sturm et al. 2006) and the submil- 
limeter galaxies (SMGs) (e.g. Alexander et al. 2005, Pope et 
al. 2008b) which are well-known sites of strong star-formation 
activity. In particular, these EDOGs could belong to the minority 
fraction of SMGs with optical-IR emission dominated by a 
dust-obscured quasar (e.g. Vignali et al. 2009) with MIR colors 
overlapping those of MIR power-law continuum AGNs (Yun 
et al. 2008; Pope et al. 2008a). Thus, even if the position 
of a source in Fig. 6a strongly depends on the criterion of 
selection adopted, EDOGs may therefore represent a mixed-bag 
population that approximately come in two flavors, i.e. "X-ray 
loud" QS02-like and SMGs-like, which essentially indicate 
different stages of quasar evolution. However, caution must 
be used before drawing any firm conclusion from these plots 
because of the uncertainties on the intrinsic Nh and the AGN 
luminosity due to the low quality of X-ray data and the use of 
photometric redshifts for most of the sources considered here. 

Fig. 6b shows the observed 2-10 keV luminosity (^?-io) as 
a function of AL^^ m for the same sources plotted in Fig. 6a. 
The solid line represents the relationship obtained from the X- 
ray sample (i.e. Eq. 2) which can be used to estimate the in- 
trinsic hard X-ray luminosity. The dotted and dashed line rep- 
resent the observed L^^-ALs^m ratio derived from Eq. (2) in 
case of a column density of Nh=10 23 and Nh=10 24 cirT 2 , re- 
spectively. Downward-pointing arrows in the plot indicate the 
3cr upper limit on L?_* for the X-ray undetected sources in our 
sample lying in the AGN region of the IRAC color-color dia- 
gram (Fig. lb), i.e. sources showing a power law spectrum and 
very red MIR colors (log(F 5 . 8 /F 3 . 6 ) > 0.2 and log(F 8 /F 4 . 5 ) > 0.2). 
Hereafter we refer to this group of 15 sources as the Red Sample. 

For these 15 X-ray undetected AGNs, we converted the up- 
per limit on the 2-10 keV flux (assuming a photon index of 
r = 1.9) into the 2-10 keV observed luminosity L?^* The in- 
trinsic L2-10 for these sources was determined from their AL^^m, 



using the relationship in Eq. (2). By the comparison between 
the observed and the intrinsic values of hard X-ray luminosity 
we were able to provide a basic estimate of the column den- 
sity distribution in the Red Sample. We found that: (i) 93% of 
these sources show a Nh > 10 22 cm -2 ; (ii) 27% of them show 
a column density Nh> 10 24 cirT 2 and are Compton-Thick AGN 
candidates; and (iii) QS02s-like objects account for 27% of the 
total (see Table 5). 

Since 75% of the Red Sample sources have X-ray effective 
exposures of ~ 5 ks, the 3cr upper limits are weakly constrained. 
Accordingly, the real values of the column density may likely 
be significantly higher than the Nh values derived here, which 
should be considered as lower limits. 

A handful of objects in the Red Sample, namely Nos. 2, 3 
and 38, were not detected in X-ray medium/deep (i.e. > 15 ks) 
observations (big downward-pointing arrows in Fig. 6b). This 
fact makes them very interesting objects since they may likely be 
more obscured than those found in the X-ray sample. In the MIR 
color-color diagram presented in Fig. lb, these EDOGs indeed 
occupy the area with the extreme red colors, i.e. log(F5 .8/P3.6) 
> 0.5 and log(F 8 /F 4 .5) > 0.5, that Polletta et al. (2008a) con- 
sider as the typical one of highly luminous and obscured quasars 
at z ~ 1. In particular, we derived an L2-10 ^ 10 44 erg s for 
sources Nos. 2 and 3 and L2-10 ~ 3 x 10 43 erg s -1 for the No. 
38, and the estimated Nh value is >10 24 cirT 2 for all these three 
sources, making them very likely Compton-Thick QS02 candi- 
dates. 

Merging the results for the Red Sample with those previ- 
ously found for the X-ray sample, we obtained a total sample 
of 38 sources for which has been possible to infer information 
on Nh and hard X-ray luminosity (see Table 5). This allows us 
to improve the statistics and get sounder insights on the X-ray 
properties^ of EDOGs, which basically turn out to host heav- 
ily obscured AGNs (~20% are consistent to be Compton-thick) 
with high X-ray luminosity ( > 40% are QS02s). 

6. Summary and Conclusions 

We exploited the multiwavelength coverage of the large area 
SWIRE survey to collect low-surface density/high luminos- 
ity (i.e. 10 12 t Lsoi ~ 10 14 L Q ) objects on the basis of 
their extreme MIR/Optical flux ratio (F2itmJpR >2000) and 
bright MIR flux density, i.e. F^pm >1.3 mJy (referred to here 
as EDOGs). Cross-correlating with publicly available XMM- 
Newton/ Chandra observations of 5 SWIRE fields we defined a 
final sample of 44 EDOGs with X-ray coverage (see Table 1 
and Sect. 2). Our study has been focused on the understanding 
of largely unexplored X-ray spectral properties of this recently- 
discovered class of sources. The main results of this paper can 
be summarized as follows: 

- All but one sources fall within the AGN region of the MIR 
color-color diagram proposed by Lacy et al. (2004). In 
particular, the vast majority of EDOGs lies in the obscured 
luminous QSO locus (Polletta et al. 2008a) characterized by 
very red colors (Fig. lb). 

- The very bulk of the redshifts, regardless of they are spec- 
troscopic or photometric, are in the range 0.7 i z ~ 2.5 
(e.g. Sect. 2.3). 



9 For the X-ray sample we considered the percentages derived using 
the lowest value in 90% confidence level for the Nh parameter obtained 
for each source. 
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- Thanks to the large area sampled (~6 deg 2 ), our EDOGs are 
on average more luminous in the MIR band by a factor of 

i 0.5-2 orders of magnitude than their equivalents in the 
CDFS and COSMOS surveys (see Fig. 3b). 

- Twenty-three sources (the so-called X-ray sample, e.g. Sect. 
4) have a positive detection either in the soft or in the hard 
band. The non-detection of the remaining sources in the 
sample can be mainly attributed to the inadequate X-ray 
coverage (most of the undetected sources indeed have an 
exposure < 5 ks). All the X-ray detected sources have a 2-10 
keV luminosity Li-\o> 10 43 erg s _1 , whereby they fall well 
within the AGN X-ray luminosity range (Fig. 4b). 

- Our analysis reveals that 50% of the sources show a column 
density strictly larger than 10 22 cm" 2 . Furthermore ~ 50% 
of them can be classified as QS02s on the basis of their 
absorption properties and X-ray luminosity (i.e. L2-io>10 44 
ergs" 1 ). 

- At least three EDOGs at z >\A were found to be very likely 
Compton-thick QS02 candidates (e.g. Sect 5.2). Moreover 
such deeply buried AGNs may be present in a sizable frac- 
tion (~ 30%, e.g. Table 5) of the X-ray undetected sources 
with extreme red colors (log(F5.g/F3.6) > 0.2 and log(Fs/F4. 5 ) 
> 0.2), as suggested by their large MIR luminosity (e.g. Sect 
5.2). They therefore deserve further investigation by future 
properly-planned, deeper X-ray observations. 

Our work represents a further step forward in completing the 
census of SMBHs through cosmic epochs. The application of our 
selection criteria to a large area survey has allowed to efficiently 
pick up a large sample of high-z sources, occupying the bright 
end of the X-ray obscured AGN population. Thus, it is a useful 
tool to complement the information derived by blind searches for 
these elusive objects performed in the pencil-beam, deep X-ray 
surveys (which mostly sample the low-redshift (z ~ 1), Seyfert- 
like AGN population) making possible the first reliable estimate 
of AGN "bolometric" luminosity function at high redshifts. 

It seems probable that EDOGs might be an early dust- 
embedded phase of the AGN activity, possibly associated with 
merger events (e.g. Georgantopoulos et al. 2008, Brodwin et al. 
2008). Accordingly, an interesting aspect that should be exam- 
ined in the next future is the FIR/sub-mm emission of these in- 
triguing sources in order to evaluate the contribution of the star- 
formation activity to their large bolometric luminosity, i.e. using 
an approach similarly to what has been done by Pope et al. 2008a 
for a sample of DOGs at low F24 /lm levels «F24^ m ) ~ 175 yuJy). 

Acknowledgements. The authors are very grateful to the referee for the exhaus- 
tive list of helpful comments and her/his prompt response. We thank for sup- 
port the Italian Space Agency (contracts ASI-INAF 1/023/05/0, ASI 1/088/06/0 
and PRIN MIUR 2006-02-5203). The authors kindly thank Simonetta Puccetti, 
Mari Polletta and Roberto Maiolino for many useful comments. Based on ob- 
servations obtained with XMM-Newton, an ESA science mission with instru- 
ments and contributions directly funded by ESA Member States and NASA. 
This research has made use of software provided by the Chandra X-ray Center 
(CXC) in the application package CIAO. This research has made use of 
the NASA/IPA Extragalactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under contract with 
the National Aeronautics and Space Administration. 

Appendix A: Notes on individual sources 

- Source No. 6. This is a narrow-line quasar at z=0.914, as 
shown in on sect. 2.3.1 on the basis of a Subaru/MOIRCS 



spectrum. This source is the brightest in our sample, with 
a 2-10 keV flux of F 2 -io~ 10~ 13 erg cirT 2 s _1 . We fitted the 
EPIC data with a continuum model consisting of a power law 
modified by an intrinsic neutral absorption component, plus 
a power law with the same slope and Galactic N H at soft en- 
ergies ( i 2 keV). The latter accounts for the "soft-excess" 
component, i.e. a mixture of photoionized/scattered emission 
superimposed upon the primary obscured continuum typi- 
cally observed in most X-ray obscured AGNs (Turner et al. 
1997; Piconcelli et al. 2007b; Cappi et al. 2006). Given the 
redshift of this EDOG, we also included in the best-fit model 
a Compton reflection component with R = 1 (where R is the 
solid angle in units of 27r subtended by the reflecting ma- 
terial) as commonly detected in the rest-frame 10-50 keV 
spectra of AGN (e.g. Risaliti 2002; Reeves et al. 2006). We 
obtained r=1.70+0.33 and N H ~ 10.8 x 10 22 crrT 2 , implying 
a 2-10 keV luminosity of L 2 -io = 5.1 x 10 44 erg s _1 . The 
most remarkable result of our analysis is, however, the dis- 
covery of an Fe Ka emission line in the X-ray spectrum of 
this QS02 (Fig. 7a). This line (E K «=6.29+0.12 keV) is re- 
quired by the data at >99.9% confidence level and shows an 
equivalent width of EW = 590+265 eV. It represents one of 
the best examples of Fe Ka line in a high-z QS02s detected 
so far. 

- Sources Nos. 12, 32 and 33. These sources (with z= 1.443, 
2.54 and 2.3 1 respectively) show most of their X-ray counts 
at E > 2 keV (e.g. Figs. 7b, 8a and 8b) and they are therefore 
promising deeply buried (possibly Compton-thick) QS02 
candidates. In particular, source No. 12 was detected with 
-40(30) net counts in the 0.5-10(2-10 keV) keV band, i.e. 
at the limit of Cash statistic reliability for the spectral fit- 
ting procedure (Sect. 4.1). Fitting these data with a simple 
spectral model consisting of a power-law with fixed photon 
index (r = 1.9) modified by an intrinsic absorption com- 
ponent, we obtained a lower limit on the column density of 
Nh > 1.4 x 10 24 cm~ 2 . The corresponding X-ray absorption- 
corrected luminosity is L2-io~ 1.4 x 10 45 erg s _1 , making 
this EDOG one of the best Compton-thick QS02 candidates 
in our sample. 

Sources Nos. 32 and 33 have only ~ 20 net counts, even if 
they were detected in deep (~70 ks) Chandra exposures, and 
almost all of them fall in the hard X-ray band. The HR of 
these sources are therefore very hard (~ 0.5, e.g. Table 3) 
and suggest that they are also highly obscured quasars. We 
derived column density values of Nh ~ 3 x 10 23 crrT 2 for 
both sources (but keep in mind the large uncertainties on the 
HR value), confirming thus the results reported in Weedman 
et al. (2006, sources A4 and A8 in their paper). Beyond pro- 
viding the spectroscopic redshifts, these authors carried out a 
detailed multi wavelength study for both sources that allowed 
to classify them as AGN-powered ULIRGs. 

- Source No. 27. This source was undetected in a ~70 ks 
Chandra exposure, being thus the only EDOG without use- 
ful X-ray information in such a deep pointing and, possibly, 
the weakest X-ray sources in our sample. Its IRAC colors 
suggest that a sizable part of the IR emission in this source 
can be powered by star formation (e.g. Sect. 2.2 and Yun 
et al. 2008). According to Rowan-Robinson et al. (2008), 
the MIR Spectral Energy Distribution (SED) is indeed dom- 
inated by the starburst component. This hypothesis is also 
supported by the large F24/F8 ratio (~20) measured for this 
EDOG. Such a value is similar to those reported for the in- 
frared luminous galaxies at z ~2 in Yan et al. 2007, for which 
both deeply embedded starburst and AGN are believed to 
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Fig. 7. a) Left panel: XMM-Newton spectrum of source No. 6 (i.e. SWIRE2 J021749.00-052306.9). The background subtracted 
spectra of EPIC PN (black), MOS1 (red) and MOS2 (green) are shown./?) Right panel: XMM-Newton spectrum of source No. 12 
(i.e. SWIRE2 J022003. 95-045220.4). The source spectrum (black stars) is superimposed to the background spectrum (red empty 
squares). The best fit models for the source (black solid line) and the background (red dotted line) are also shown. 
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Fig. S.a) Left panel: Chandra spectrum of source No. 32 (i.e. SWIRE2 J104409. 95+585224. 7). b) Right panel: Chandra spectrum 
of source No. 33 (i.e. SWIRE2 J104528. 29+591326. 6). Black stars (red empty squares) represent the source (background) binned 
counts, with 2 (15) counts per bin. 



significantly contribute to the dust heating (e.g. Polletta et 
al. 2008a). 
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